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We propose an explanation of the LSND signal via quantum-decoherence of the mass states, 
which leads to damping of the interference terms in the oscillation probabilities. The decoherence 
00 ■ parameters as well as their energy dependence are chosen in such a way that the damping affects 

only oscillations with the large (atmospheric) Am 2 and rapidly decreases with the neutrino en- 
fN) ' ergy. This allows us to reconcile the positive LSND signal with MiniBooNE and other null-result 

I/") ■ experiments. The standard explanations of solar, atmospheric, KamLAND and MINOS data are 

not affected. No new particles, and in particular, no sterile neutrinos are needed. The LSND 
signal is controlled by the 1-3 mixing angle 6*13 and, depending on the degree of damping, yields 
0.0014 < sin 2 #13 < 0.034 at 3d. The scenario can be tested at upcoming #13 searches: while the 
^ comparison of near and far detector measurements at reactors should lead to a null-result a positive 

signal for $13 is expected in long-baseline accelerator experiments. The proposed decoherence may 
partially explain the results of Gallium detector calibrations and it can strongly affect supernova 
neutrino signals. 
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I. INTRODUCTION 



Reconciling the LSND signal for^i^ transitions [l| with the by now fully established 

on the one hand side and bounds on oscilla- 
13j on the other is a long-standing problem in neutrino physics. 

failed to confirm the oscillation interpretation of 



U1 "JL ~^ !f£ 11 

existence of neutrino oscillations 
tions @, 0, 0, 0, 
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Recently the MiniBooNE experiment 
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the LSND signal [15j, putting further pressure on its interpretation. Usually the LSND 
result is considered as an indication for sterile neutrino oscillations, despite difficulties of 
the corresponding models to explain the global data 16] including cosmological observations 
(see Ref. 17] for a recent analysis including the MiniBooNE results). Triggered by these 
problems many ideas and scenarios have been proposed in order to explain LSND, some of 
them involving very exotic physics. This includes sterile neutrino decay la 
of the CPT [20] and/or Lorentz 21] symmetries, mass- varying neutrinos 22 



violation 
short-cuts of 

sterile neutrinos in extra dimensions 23], a non-standard energy dependence of the sterile 
neutrino parameters j^], or sterile neutrinos interacting with a new gauge boson 25 . 



In the present paper, we revisit the possibility that the origin of the LSND signal might 
be quantum decoherence in neutrino oscillations [26j, |27]. Such effects can be induced by 
interactions with a stochastic environment; a possible source for this kind of effect might be 
quantum gravity 



28, 29, 30 



see Ref. [31] for a recent discussion. We take a phenomenolog- 
ical approach and determine the form and magnitude of the new effects from observations 
without reference to possible origins of the decoherence. Neutrino oscillations, being a quan- 
tum interference effect over macroscopic distances, provide a sensitive test for decoherence. 



quantum decoherence has been explored for 
KamLAND 33 .I34] . future long-baseline ex- 



The possibility to use neutrinos as a probe o 
atmospheric neutrinos |32j , solar neutrinos [33 
periments [35j, |36J , and neutrino telescopes [3 
which have been considered to search for decoherence include, for example, KK 29_|, |39_|, |40 
and BB 4l|, |42J oscillations or neutron interferometry 43] . Throughout this paper we will 



See also Ref. [381] . Other quantum systems 



assume that quantum decoherence affects only the neutrino sector. 

Previous attempts to explain the LSND signal by quantum decoherence 26j, |27J seem to 
be in conflict with the present data. Indeed, for the set-up of NuTeV with a baseline of ~ 1 



km and the average energy of 75 GeV the model in [27J predicts P( 



u e ) = Piv^ 



v,. 



3 x 10 3 while the one in 26j yields -P(y M — > v e ) — 0.2. Both of these predictions strongly 
violate the upper bound on the neutrino and anti-neutrino oscillation probabilities from 
NuTeV: P{v^ -»■ u e ), F(p„ -»• v e ) < 5 x 10~ 4 (90% C.L.) 0- Furthermore, the model of [26 
(where in addition to decoherence, also CPT-violation is also introduced) cannot account 
for the spectral distortion in the anti- neutrino signal observed by KamLAND. The scenario 



of [27] is also disfavored by the absence of a signal in KARMEN [9], NOMAD 12] and 



MiniBooNE 
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In the present paper we propose a different set of decoherence parameters 
with a special energy dependence. As a result, the problems outlined above can be avoided. 
In our scenario, the decoherence effects become more significant with decreasing the energy; 
thus, we refer to this scenario as "soft decoherence" . 

The paper is organized as follows. In sec. [Til we describe our scenario of the decoherence 
effect and discuss how the LSND signal is explained. In sec. UTIJ we show that by selecting 
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suitable parameters, within this scenario the LSND result can be reconciled with all other 
oscillation results (both positive and negative). In sec. IIV} we discuss how this scenario 
can be checked in forthcoming and future experiments. Discussions and conclusions will be 
presented in sec. [V] 



II. THE SCENARIO: DECOHERENCE AND THE LSND RESULT 



Soft decoherence 



Let us describe the neutrino system by the density matrix p in the mass state basis, fj, i = 
1, 2, 3. The decoherence effects in the evolution of the density matrix can be parameterized 
by introducing a new term T>[p] as 



^ = -i[H,p]-V[p}. 



(1) 



This term violates the conservation of Tr(p 2 ) and, hence, leads to the evolution of pure states 
into mixed states. As mentioned in introduction, such a term can be induced by interaction 
with a stochastic environment; a possible source for this kind of effect might be quantum 
gravity 28|, |29|, |30(. The form of the operator T>[p] can be constrained by imposing some 



general requirements on the evolution of the system. First, complete positivity implies the 
so-called Lindblad form for T>[p] jiil |45| : 



V[p}=J2[{P,D n Di}-2D n pDl} 



(2) 



where D n are some general complex matrices. This form arises from "tracing away" the 
dynamics of the environment 3C|, |46) . Second, with general complex D n unitarity is violated, 
i.e., dTr(p)/dt can be nonzero. Therefore, we require that D n are Hermitian, = D n . 
In addition to dTi{p)/dt = 0, the Hermiticity of D n guarantees that the entropy [i.e., 
S(p) = — Tr(plnp)] cannot decrease 45|, |47|. Finally, we require that the average energy of 
the system, Tr(pH), is conserved. It is straightforward to check that this can be achieved by 
demanding [H, D n ] = 0. In other words, unitarity and conservation of the energy-momentum 
imply that D n and H can be simultaneously diagonalized. In the neutrino mass basis, we 
can therefore write 



H = Diag[/i 1; h 2 , h 3 ] and D n = Diag[d ni i, d n>2 , d n - 



(3) 



where hf = p 2 + m 2 , and d n ^ are real quantities of dimension [mass] 1 / 2 whose energy- 
dependence is unknown. In this paper we adopt a phenomenological approach and determine 
d n> i from observations without discussing their possible origins. 

Solving the evolution equation Eq. ([!]), with H and D n given in Q we find 



Pit) 



Pn(0) pi 2 (0)e-^ 12 - iAl2 )* pi3(0)e^ 13 - iAl3 )*' 

p 2 i(0)e-(^ A21 )' p 22 (0) p 23 (0)e-^ :i - iA ^ t 

Lp3i(0)e^- iA31 )' p 32 (0)e-^ 2 - iA32 )* p 33 (0) 



(4) 



3 



where Pij(0) are the elements of the density matrix at the initial moment, 

/\t7?,^ - 

lij = ^2(d n ,i ~ d n j) 2 and A j{ = hj - hi « . (5) 

n v 

Notice that 7^ = 7^ whereas Ay = — A^, and the diagonal elements of p do not depend on 
time. 

Let us consider the transitions between the flavor states, v a = Yli^ai^ii where U a i are 
the elements of the PMNS mixing matrix. The probability of finding a neutrino with flavor 
(3 is given by (vp\p\vp). Hence, the oscillation probability v a — > up in vacuum is equal to 

P aP = (Mp^WWp) = E u h u v p§°(*) > ( 6 ) 

ij 

where p\j(t) is given by Eq. (@J with Pij(0) = Pi" (0) = U a iU*j, which corresponds to the 
initial state v a . 

In this paper, we consider the most economic scenario that describes all the data. As we 
will show in the following, only one matrix D n with 

d x = d 2 ^ d 3 (7) 

is sufficient. Note that a similar pattern exists between neutrino masses (mi ~ 7712 7^ ^3), 
so it will be inspiring to build a model that links the patterns of H and D. Eq. © leads to 

7i 2 = and 7 = 713 = 7 32 . (8) 

In the ranges of energy and baseline (L) for which A 2 ±L = Am^I / '(2E„) <C 1 , the oscilla- 
tions due to AmJ-L can be neglected, and Eq. (jSJ) yields 



P Me ( 7 ,L) = P eAt ( 7 ,L) = 2|^ t 3| 2 |f/ e 3| 2 [l-e-^cos(A 31 , 
P ee ( 7 , L) = 1 - 2|f/ e3 | 2 (l - |t/ e3 | 2 ) [1 - e^ L cos(A 3 iL)] , (9) 
P^( 7 ,L) = 1-2|^ 3 | 2 (1-|M 2 ) [l-e- 7L cos(A 31 L) 



Notice that although the new term T>[p] explicitly breaks the time reversal symmetry, 
still within the framework of the two-neutrino oscillation the equality P eM = P Me holds. In 
principle, the decoherence effects can give rise to the CPT violation; however, in this paper 
we assume that the decoherence effects in the neutrino and antineutrino sectors are the 
same. 

The energy dependence of the decoherence parameter 7 in Eq. (jSJ) is not known; it should 
follow from a microscopic theory of decoherence. In the absence of such a theory we assume 
a power law: 7 oc E~ r , and for convenience parameterize it as 

,2 //in A/T^AT" \ r ~ 1 



p 

7 = a 



where 40 MeV is the typical neutrino energy in LSND. We will estimate the allowed ranges 
of parameters r and p 2 in sec. IHIt and in particular, show that all the data can be described 
if 

r = 4, (11) 



4 



which we will use as the reference value in our estimations. The fast decrease of 7, and 
consequently, the decoherence effect with energy, is the key feature of the proposed scenario 
which allows us to reconcile the LSND result with results of other experiments. Notice that 
for any value of r 7^ 1, the Lorentz symmetry is explicitly violated. 



B. Explaining LSND events 

Let us now discuss the interpretation of the positive LSND result through soft deco- 
herence. For the LSND parameters, E ~ 40 MeV and L ~ 30 m, the oscillation phase 
A 31 L ~ 5 ■ 10~ 3 , so that the oscillation effect is negligible and according to (Q the appear- 
ance probability is equal to 

P Me ( 7 , L) = 2|f/ M3 | 2 |f/ e3 | 2 (1 - e-' L ) « \U e3 \ 2 (1 - e^ L ) . (12) 

Hence, the LSND signal is determined by the 1-3 mixing and the degree of decoherence given 
by the factor (1 — e~ lL ). In the case of strong decoherence, 7L ^> 1, this factor converges 
to its maximum and 

P Me ( 7 ,L)^|[/ e3 | 2 . (13) 

The LSND signal is simply given by the 1-3 mixing parameter, and therefore the LSND 
probability provides a lower bound on the 1-3 mixing: 

|^ e3 | 2 > P^ SND = (2.6 ± 0.8) • 10" 3 . (14) 
For a small decoherence effect, 7L <C 1, we obtain 

P Ate ( 7 ,L)^|f/ e3 | 2 7^. (15) 

Therefore \U e ^\ 2 can be as large as the upper bound obtained from the reactor experiments 
and other observations, (|f/ e 3| 2 ) upper - I n this case the LSND result gives a lower bound on 7 
at typical LSND energies: 

pLSND 

7 - (|f/ e3 | 2 ) UPPe ^LSND 1 ] 

In Fig. [1] we show the Lj E dependence of the excess of the LSND events. The points with 
error bars are the LSND data and the histograms show the prediction of the decoherence 
scenario (added to the background) for various values of the power r, see Eq. ( flOj) . To draw 
each histogram, the corresponding best fit value for fi 2 is assumed. For comparison we show 
the spectrum expected due to oscillations in the model with a light sterile neutrino. The 
decoherence reproduces the observed spectrum quite well. It leads to a softer energy spec- 
trum than the oscillation spectrum. Furthermore, with increase of r the softness increases 
and the maximum shifts to lower energies. Unfortunately these differences are small and 
with the LSND statistics and uncertainties it is not possible to distinguish decoherence and 
oscillation effects or substantially restrict r for r > 4. For 7L ~ 1, L ~ 30 m and typical 
LSND neutrino energies E v ~ 40 MeV, the relevant scale for the parameter \i introduced in 
Eq. (TTO]) is /x 2 ~ 0.1 eV 2 . 
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FIG. 1: Spectrum of v e excess events in LSND. The histograms correspond to the prediction of the 
decoherence scenario for various values of r, see Eq. (Tl0|) . The best fit value for /x 2 in each case is assumed. 



Thus, in the proposed scenario the LSND — > u e signal is explained by the decoherence 
of the mass state z/3, whose mixing with and v e are given by cos #13 sin #23 and sin #13, 
respectively. The oscillation effect is negligible and plays no role here. Linking the LSND 
signal with the last unknown mixing angle, #13, is an exciting feature of the model. We will 
discuss the implications of this relation in sec. IIVI 



III. RECONCILING LSND WITH OTHER OBSERVATIONS 

To simultaneously accommodate the positive LSND signal and the lack of any evidence for 
flavor transitions from other short-baseline experiments, the value and energy-dependence 
of 7 have to be properly chosen. In what follows, we demonstrate that since the neutrino 
energy and baseline for each experiment are different, such a choice is possible without 
affecting the successful description of solar, atmospheric, KamLAND, and MINOS data in 
terms of neutrino oscillations. 



A. Short baseline and reactor experiments 

In short-baseline experiments, Am^L/(2E u ) <C 1 and we can neglect oscillations. There- 
fore, for a given power r there are three parameters to describe these experiments: /i 2 , 
\U e3 \ 2 = sin 2 13 and \U^ 3 \ 2 = cos 2 6*13 sin 2 #23 [see Eq. fl9])]. Hence, this scenario in- 
volves the same number of parameters as the (3+1) sterile neutrino oscillation schemes. 
We have performed a fit to the data from LSND [l[ (decay-at-rest data), KARMEN Q], 
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FIG. 2: Allowed regions in the plane of |^ e 3| 2 and the decoherence parameter for r = 4. Shaded regions: 
global data at 90% and 99% CL. Curves: 99% CL regions from LSND, KARMEN, MiniBooNE and reactor 
data (Bugey, Chooz, Palo Verde). The star marks the global best fit point. We marginalize over sin 2 #23 
taking into account the constraint from atmospheric neutrinos. 



MiniBooNE 0, H] , NOMAD Q, CDHS (h}, and the reactor experiments Bugey [n 



17, 49 . Results of the fit are 



Chooz jjj, and Palo Verde |8j. For technical details, see Refs. 
shown in the ^ 2 -\U e ^\ 2 plane in Fig. [2j 

The LSND signal can be reconciled with the null-result from KARMEN (which has 
the same neutrino energy as LSND) due to the somewhat shorter baseline in KARMEN 
(L ~ 18 m) and the exponential dependence of the decoherence effect on the distance. The 
solid and dashed curves in Fig. [2] show the allowed regions from LSND and KARMEN, 
respectively. At 99% C.L., there is a significant overlap left, though there remains some 
tension between the two results. Notice that a similar tension exists in other scenarios 



developed to explain the LSND results such as the case of oscillations [Ha, [50] or sterile neu- 
trino decay 19( , the latter having the same exponential L dependence as in the decoherence 
model. 

The baseline and energy of neutrinos at MiniBooNE are one order of magnitude larger 
than in LSND. As a result, the oscillation phase Am 2 L/(2E l/ ) for the two experiments are 
of the same order and oscillations can be neglected. We find that for r > 2, at MiniBooNE 
7L <C 1 and therefore the decoherence effects are negligible, rendering the new effects at 
MiniBooNE unobservable and explaining the non-observation of an appearance signal. As 
visible in Fig. [21 for MiniBooNE the decoherence becomes important only for yU 2 > 20 eV 2 , 
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and is completely negligible for values of fi relevant for LSND. Similarly to the case of 
standard oscillations, the soft decoherence cannot account for the event excess observed in 



MiniBooNE below 475 MeV. Therefore, following Ref. [14j, we rely on a yet to be identified 
explanation of this excess and use only the data above 475 MeV in the analysis. 

At the reference value r = 4 any decoherence effect for short baseline experiments with 
energies substantially larger than 40 MeV is suppressed. This implies no flavor transitions in 
the "high energy" experiments {E v > 0.5 GeV) MiniBooNE, NOMAD, NuTeV, and CDHS. 1 
Therefore, our scenario is trivially consistent with the null-result of these experiments. 

The energies of reactor neutrinos are relatively low (of order 4 MeV). As a result, for 
these neutrinos 7 is quite large: 

or fi 2 /4MeV\ 4 l 
7~2.5 — - — ^ — - — cm 1 . (17) 
[ 0.2 eV 2 \ E v ) y ! 

Thus, already after a few centimeters coherence is completely lost and the survival proba- 
bility in Eq. (Q becomes independent of the energy and baseline: 

pir tor ~i-2|f/ e3 | 2 (i-if/ e3 i 2 ). (is) 

Thus, the Chooz bound on P ss , which is derived by comparing the measured v e flux at a 
distance of 1 km from the source with the estimated flux from the consideration of the power 
of the reactor can be directly translated into an upper bound on |t/ e3 | 2 . The upper bound 
on |£/e3| 2 from the combination of Chooz, Bugey and Palo Verde data is shown in Fig. [21 
We find \\J e3 \ 2 < 0.04 at 3a (1 d.o.f.). 

The shaded regions in Fig. [2] show the results of the global analysis of the short-baseline 
experiments. We include sin 2 #23 as a free parameter in the fit, taking into account the 
standard constraint from the Super-Kamiokande atmospheric neutrino data {2J, which is 
hardly affected by decoherence (see next section for a more detailed discussion). For r = 4 
we find the best fit point at 



\U e3 \ 2 = 6.1 x 10~ 3 , sin 2 6 23 = 0.5, /i 2 = 0.27eV 2 , (19) 

with Xmin = 89/ (107 — 3) d.o.f. For these values of the parameters, the averaged probability 
in LSND equals -P^ e = 2.3 • 10~ 3 which is within one sigma of the experimental value 
(2.6 ± 0.8) • 10~ 3 . The allowed range for \U e3 \ 2 is 

(1.4)2.3 x HT 3 < \U e3 \ 2 < 2.1 (3.4) x 10~ 2 at 2 (3)a (1 d.o.f.) . (20) 

Here the lower bound follows from the LSND result in the limit of strong decoherence, see 
Eq. ffl4j) . The upper bound comes from the reactor experiments. 

Let us evaluate the quality of the fit in more detail, and compare it to the case of sterile 



neutrino oscillations in a (3+2) scheme [5jJ. To this aim we divide the data into sub-sets 



1 A similar strategy to reconcile low and high energy short-baseline data has been pointed out in Ref. 24 1 
in the context of sterile neutrinos with a non-standard energy dependence. 
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(3+2) oscillations decoherence 

Data sets xp G / d -°- f - PG Xp G / d -°- f - PG 

LSND vs NEV 21.2/5 0.08% 8.6/2 1.4% 

App vs Disapp 17.2/4 0.18% 0.6/2 74% 

TABLE I: Consistency of LSND versus all other no-evidence short-baseline data (NEV), and appearance 
versus disappearance short-baseline data for (3+2) oscillations and the decoherence model. We give Xpg 
according to Eq. (f2"Tj) and the corresponding probability ("PG"). The results for (3+2) are taken from Tab. 3 
of Ref. [171. 

and check the consistency of these data sets by using the so-called Parameter Goodness-of-fit 



(PG) criterion 1_6|, [52j. It is based on the x 2 function 

XPG Xtot,min ,min ) 

(21) 

i 

where Xtot min ^ s the x 2 minimum of all data sets combined and xf jmin is the minimum of the 
data set i. This x 2 function measures the "price" one has to pay by the combination of the 
data sets compared to fitting them independently. It should be evaluated for the number of 



d.o.f. corresponding to the number of parameters in common to the data sets, see Ref. 52 
for a precise definition. 

First we test the consistency of LSND with all the other null-result short-baseline ex- 
periments (NEV). The numbers given in Tab. [I] show that for (3+2) oscillations LSND is 
consistent with NEV only with a probability of 0.08%, whereas in the decoherence scenario 
the probability improves to 1.4%. Note that in the (3+2) case short-baseline data depend 
on 7 parameters, whereas for decoherence only 3 parameters are available to fit the data. 2 
The reason for the still relatively low probability of 1.4% is the aforementioned tension 
between the LSND and KARMEN results. In order to illustrate this effect we perform a 
second test, dividing the data into appearance experiments (LSND, KARMEN, NOMAD, 
MiniBooNE) and disappearance experiments (CDHS, Bugey, Chooz, Palo Verde). In this 
approach, LSND and KARMEN data are summed into the same data set and therefore, 
by assumption they are taken to be consistent. As a result, the remaining tension between 
them does not show up in the PG test. Note that it is reasonable to combine LSND and 
MiniBooNE into the same data set, because in both cases considered here they are consis- 



tent: for (3+2) oscillations they can be reconciled [17( by invoking CP violation [19|, 153 
whereas in the decoherence scenario the energy dependence of 7 guarantees the null- result 
of MiniBooNE. As shown in Tab. IH we find an excellent fit in the decoherence model (PG 
of 74%), whereas (3+2) oscillations suffer from a sever tension between appearance and 
disappearance experiments, allowing for compatibility with a probability of only 0.18%. 

In summary, the soft decoherence proposed here provides an excellent fit to short-baseline 
experiments, allowing for full consistency of LSND and MiniBooNE, as well as appearance 



2 As before we work at fixed r = 4 and do not consider the energy exponent as a free parameter. 



9 



and disappearance experiments. Only the well-known tension between LSND and KARMEN 
remains unresolved. 



B. Other phenomenological consequences 

In this section, we show that our scenario is compatible with the standard description 
of the solar [5], KamLAND Q, atmospheric [2], K2K [4] and MINOS 0] data in terms of 
neutrino oscillations. 



Solar and long-baseline reactor neutrino data: In Ref. 33J, decoherence effects 
on solar and KamLAND neutrino data have been studied (see also Ref. 34 for the case 



of KamLAND only). The decoherence scenario in 33j differs from the one in the present 
paper: While we take 712 = (d\ — d 2 ) and discuss the effects of 7 = {d\ — d 3 ) 2 , the 



authors of |33| focused on the effects of 712 = {d\ — d 2 ) 2 . Our assumption d\ = d 2 in Eq. ([7]) 



ensures that oscillations in the 1-2 sector are not affected by decoherence. This leaves 
the dominant oscillations due to Am^ and 812 unchanged and guarantees the standard 
oscillation explanations for the solar and KamLAND data. Within our scenario, the effects 
of the damping factor on the solar neutrino flux and KamLAND neutrinos are suppressed 
by I t/ e 3 1 2 . The current uncertainties do not allow to resolve the effects of |t/ e 3| 2 . Moreover, 
in these experiments oscillations due to Am^ are completely averaged out and therefore 
decoherence effects in the 1-3 sector are unobservable. 

We can use results of | 3 31] to put an upper bound on 712 = {d\ — d 2 ) 2 . Writing 712 = 
7o(l GeV I E u ) r , in Ref. [331] bounds on 70 have been derived from solar and KamLAND data 
assuming r = 0, ±1,±2. Extrapolating these results to r = 4, one finds 70 < 10~ 32 GeV. 
Using the parametrization shown in Eq. (flu]) , one has 70 = 6.4 • 10~ 23 GeV (/^/eV 2 ), and 



fx 2 12 < 10- 10 eV 2 « /i 2 



(22) 



Atmospheric neutrinos: First, we note that because of the smallness of |£7e3| 2 , the 
decoherence effects do not considerably change the v e flux at low energies where the Earth 
matter effect can be neglected. For high energies (multi-GeV sample) where the 1-3 mix- 
ing is enhanced the decoherence effect becomes negligible. However, the muon neutrino 
disappearance probability can be significantly affected. In order to have a detectable ef- 
fect of decoherence, the neutrino energies and baselines should be in a range for which 
7L > 1 and A 31 L < 2n. For A 31 L ^> 2tt, regardless of the value of 7, the interference 
term e~ lL cos(A 31 L) averages out and the sensitivity to the new effects is lost. Putting the 
two conditions together, we find that the effects can be noticeable only for E v < 400 MeV. 
On the other hand, for E u < 200 MeV, the produced muons cannot be detected in Super- 
Kamiokande. 

In Fig. [3] we show dependence of the survival probability — > z/ M on the neutrino zenith 
angle for different energy ranges. In the lowest energy range, 200 MeV < E u < 400 MeV, 
the oscillation length L osc = 2TrE u /Am 2 l ~ 100 km. Hence, the condition A31L ~ 7r implies 
that only for neutrinos arriving from above (for zenith angles smaller than 90°) the damping 
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FIG. 3: The zenith angle dependence of the survival probability relevant for atmospheric neutrinos 
in three different energy intervals for r = 3,4 and for the case of oscillations without decoherence. The 
horizontal dashed lines show the corresponding probabilities averaged over all zenith angles. 



effects are significant. As follows from Fig. [31 in the vertical direction for r = 4 and r = 3 the 
effect can reach 10% and 30%, correspondingly. However, for these low energies the direction 
of the initial neutrino is not related to the muon direction and hence, the distribution of 
/i-like events is averaged over the zenith angle. As a result, the decoherence leads to a 
zenith angle independent decrease of the sub-GeV /i-like events. The size of the effect is 
illustrated by the horizontal dashed lines in Fig. [3], which correspond to the zenith angle 
averaged survival probability. For 200 MeV < E v < 400 MeV, we find a suppression of 
5% (18%) for r = 4 (r = 3). This effect can be detected as a decrease of the ratio of the 
/i-like to e-like events. We find that for r > 4 this decrease is below the 5% experimental 
uncertainty on this ratio 0]. Essentially this effect determines the lower bound on r. Let us 
note that at these low energies, there is some excess of e-like events in Super-K, which can 
be interpreted as renormalization and deficit of /i-like events and the latter can be explained 
by the decoherence in our scenario. 

In the energy interval 400 MeV < E v < 1.3 GeV the event suppression due to the 
decoherence effect is below 5% even for r = 3. The decoherence leads to a flattening of the 
zenith angle dependence. However, the averaging over the zenith angle in the /i-like events — 
though incomplete — is still strong. So uncertainties in the extraction of the zenith angle in 
the Super-Kamiokande experiment do not allow us to identify the effect of decoherence for 
r > 3. In the multi-GeV range the decoherence effect is strongly suppressed. In the Soudan 
experiment 54( , in addition to the energy of the produced muon, the recoil energy of proton 
(in quasi-elastic interaction) is also measured. As a result, the zenith angle of the incoming 
neutrino can be deduced. However, in this experiment, the statistical error is larger than 
10% and the decoherence effects cannot be resolved. We conclude that for r > 4 our scenario 
is consistent with the oscillation interpretation of atmospheric neutrino data. 3 



We thank Michele Maltoni for communication on this point. 
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Long-baseline accelerator experiments: At MINOS 0] experiment with E u >2 GeV 
and L = 730 km, the damping effect is very small, (1 — e~ 7i ) ~ 7L < 0.03, and the present 
uncertainties do not allow to resolve the effect. In the K2K experiment [4( the energies of 
the neutrinos are smaller and the decoherence factor is slightly larger. However, still the 
statistics below GeV is too low for the experiment to be sensitive to the soft decoherence. 

In Ref. 32], the decoherence effects in the atmospheric and K2K neutrino data have been 
studied. However, the bounds derived in [32J do not apply here because we assume a steeper 
decrease of 7 with energy. 

Radioactive source experiments: In calibrations of the Gallium solar neutrino de- 
tectors SAGE 55| and GALLEX/GNO 56j with artificial radioactive 51 Cr and 37 Ar sources 
deficits of the signals have been reported in Ref. [55| : the weighted average of the ratio of 



the observed to expected event numbers equals i?Ga = 0.88 ± 0.05 [55j. This result was in- 



terpreted as an indication of electron neutrino disappearance due to oscillations into sterile 
neutrinos 57[. For the energies of the radioactive sources employed in these experiments, 
E v ~ 0.8 MeV, the damping parameter 7 is huge [see Eq. ffTTj) ] and the decoherence length 
~ I/7 ~ 5 ■ 10~ 4 cm. So, complete decoherence occurs over a small fraction of a millimeter 
and the survival probability is given by Eq. ([TBI . For the best fit value of |t/ e 3| 2 shown in 
Eq. (TT9l . we find a small effect: P ee ~ 0.99. However, if \U e z\ 2 is at its 3a upper bound given 
in Eq. ([201 we obtain P ee ~ 0.93. The latter is within the la range of Rc a and therefore the 



results of the calibration experiments can be at least partially explained. 



IV. FUTURE TESTS OF THE SCENARIO 
A. Reactors versus accelerators 

The possibility to test decoherence effects at future long-baseline experiments has been 



discussed in [35|, |36|. In our scenario, an explanation of the LSND signal implies a lower 
bound on the mixing angle #13. Therefore, upcoming oscillation experiments aiming at the 
measurement of this mixing angle will provide a crucial test. 



The next generation of reactor experiments like Double-Chooz 58j, Daya Bay |59j and 



Reno [60] will search for # 13 by comparing the anti-neutrino flux measured at a near and far 
detector. As discussed in sec. IIII A\ a remarkable consequence of our decoherence scenario 
is that for baselines larger than a few centimeters the interference disappears, and P ss does 
not vary with L or E u , see Eq. ffTBl . As a result, the comparison of signals in near and 
far detectors at reactors will not reveal oscillations. The model predicts a v e flux reduction 
already at the near detector with respect to the initial flux emitted from the reactor. How- 
ever, establishing this reduction would rely on the ability to determine the original flux with 
better than 1% accuracy. This seems difficult to achieve. 



In contrast to the reactors, in the future long-baseline accelerator experiments T2K [61 
and NOi/A [621 the damping effect will be quite small, (1 — e _7L ) ~ 7L ~ 0.02, due to 
smallness of 7 in the GeV energy range. Consequently, the sensitivity of these experiments 

fie 



to \U e3 \ 2 through measurements of P„ e will be basically unaffected by decoherence. The 
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sensitivity of T2K and NOi/A is in the range sin 2 2#i 3 > 0.01, and therefore these experiments 
can probe the best fit point of our scenario, Eq. (fl9l) . and large part of the allowed range 
Eq. (I20l . though the lower end of the 3a interval for |£7 e 3| 2 may escape detection. Hence, 
by comparing the results of T2K and the upcoming reactor experiments, the validity of the 
present scenario can be tested. The following three situations are of particular interest. 

• Comparing the flux at far and near detectors, the reactor experiments such as Double- 
Chooz and Daya Bay would establish a non-zero value of \U e 3\ 2 m agreement with the 
value extracted from T2K and NOi^A measurements of P^ e . In this case, our scenario 
would be ruled out because, as discussed above, P g g at the near and far detectors 
should be the same. 

• Neither the reactor experiments nor the long baseline experiments find any evidence 
for nonzero \U e3 \ 2 and put an upper bound of 0.0025 on its value (sin 2 26*i 3 < 0.01). In 
this case, the constraints on our scenario become quite tight, and the allowed region 
for |£/ e3 | 2 shifts to the lower end of the 3cr interval given in Eq. (!20j) . This shifts yU 2 
to larger values, of order 1 eV 2 (see Fig. [2]) and increases the tension between LSND 
and KARMEN data. In order to fully rule out the model, the 3a bound on |?7 e3 | 2 has 
to be pushed below 0.0014 (sin 2 26* 13 < 0.0056), which probably requires to go beyond 
the initial phases of T2K and NOz/A. 

• While comparing the reactor neutrino fluxes at near and far detectors reveals no ev- 
idence for missing z/ e , T2K and/or NOz/A consistently report a value of \U e ^\ 2 in the 
range (0.0014,0.034). This situation cannot happen within the standard oscillation 
scenario. Thus, such an outcome can be considered as a strong hint in favor of our 
scenario. A way to confirm or refute this hint is to compare the measured reactor 
neutrino flux with the original flux estimated from the power considerations, as it had 
been done in the analysis of the CHOOZ data. Of course, to do this the systematical 
uncertainties in the flux estimations have to be overcome. 

Currently the MiniBooNE experiment is taking data in the anti-neutrino mode. Since 
our model invokes neither CP nor CPT violation, the prediction for anti-neutrinos is the 
same as for neutrinos. Hence one expects a null-result for MiniBooNE anti-neutrino run. 



For the low energy experiment proposed in Ref. [63j using the LENS detector, the situation 
is similar to the one in the Gallium calibration experiments mentioned above. After very 
short distances decoherence sets in, leading to a constant event suppression according to 
Eq. ( fl8l) . No distance dependent effect would be observed, and the measurement has to rely 
on the comparison of expected and predicted numbers of events, which might be difficult 
due to the normalization uncertainties. 

The phase I of T2K is planned to be followed by a phase II which can probe the effects 
of Airijj and, for relatively large values of \U e ^\ as predicted in our model, measure the 



Dirac CP- violating phase 6jj . The energy of the neutrino flux in the second phase of T2K 
will be around 750 MeV for which 7 ~ Am^/^, Thus, the decoherence parameter is 
(1 — e~ lL ) ~ 7L ~ Amjj/Am^ ~ 0.03. Hence, one expects a distortion of the energy 
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spectrum at the level of a few percent. The effect might be difficult to observe in the 
appearance signal, since in this case, the number of events expected in T2K-II is of order 
of 1000, i.e., the statistical error is a few percent. Moreover, from Eqs. (jHJ), we observe 
that the decoherence effects on the appearance probability, P^ e , is further suppressed with 
a factor of |t/ e 3| 2 . However, a spectral distortion due to decoherence is also expected for the 
Up disappearance signal, where statistics is much larger. 

Future projects for superbeam (SPL) or Beta Beam experiments from CERN to a mega- 



ton scale detector in Frejus |64j have very good sensitivity to probe our scenario. These 
experiments will have neutrino energies roughly a factor two smaller than T2K, which en- 
hances 7 by a factor 2 r and therefore, jL ~ 0(1). Thus, the experiments would operate 
in the regime that the decoherence effect is significant and this would lead to very different 
spectral signatures as compared to standard oscillations. 



Ref. [65] suggests to use so-called Mossbauer neutrinos to measure 6*i 3 . The energy of 
such neutrinos is low (for example, for neutrinos from Tritium decay, E u = 18.60 keV). For 
such neutrinos, our scenario predicts Eq. (fl8l) . while within the standard oscillation scenario, 
P ee = l- 4|[/ e3 | 2 (l - \U e3 \ 2 ) sin 2 A 3 iL/2 [66(]. Considering that the energy spectrum of the 
Mossbauer neutrinos is monochromatic, it will not be possible to discriminate between the 
two scenarios by studying the energy dependence of P ee . However, by measuring the flux at 
several distances it will be possible to make a distinction. If the soft decoherence is realized 
in nature and the measurement is done at a distance L with sin 2 A31L/2 7^ 1/2, neglecting 
the decoherence effects will cause a disagreement between the results of T2K and NOz/A 
with this measurement. 



B. Decoherence in matter; Supernova neutrinos 

In matter, the total Hamiltonian includes the interaction term described by the matrix 
of potentials V: H — > H + V. In the neutrino mass basis V is non-diagonal, and therefore 
[H, D n ] 7^ 0. The fact that the decoherence matrix does not commute with V and conse- 
quently with the total Hamiltonian leads to a new effect: statistical equilibration of flavors 
and masses. This means that after sufficient time (t ~ ([V, D])^ 1 ^ 2 ), p converges to unit 
matrix times a normalization factor and consequently, in the case of two neutrino mixing, 
the probabilities of finding neutrinos with the mass mi and ni2 in the course of evolution 
converge to Pi = P2 = 1/2. Similarly for mixed flavors as a result of long enough evolution 
P e = P x = 1/2, where x is some combination of and v T neutrinos with which v e mixes. 
This phenomenon is similar to the flavor equilibration in the presence of mixing and inelastic 
collisions which destroys coherence. Such a type of equilibration has been considered for the 
active-sterile neutrino oscillations in the Early Universe. 

In our scenario non-trivial interplay of the decoherence and matter effect should take 
place for solar, atmospheric and supernova neutrinos. No significant effect is expected for 
the solar neutrinos inside the Sun as well the Earth and also for the atmospheric neutrinos 
inside the Earth. Indeed, for solar neutrinos due to low energies the matter effect on 1-3 
mixing is negligible. For atmospheric neutrinos inside the earth the matter effect on 1-3 
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mixing is substantial in the multi GeV range where decoherence is negligible. In contrast, 
for supernova neutrinos (E = 5 — 40 MeV, and huge densities) both matter and decoherence 
effects are strong. 

Let us estimate qualitatively the decoherence effect on supernova neutrinos. In our sce- 
nario the 1-3 mixing should be relatively large so that without decoherence the conversion in 
the H-resonance region (due to 1-3 mixing and mass splitting) should be highly adiabatic. In 
the case of normal mass hierarchy that would lead to the transition v e —>■ z/ 3 inside the star. 
As a consequence, no earth matter effect is expected in the neutrino channel. In contrast, in 
the presence of decoherence and equilibration of mass states only half of z/ e 's will end up as 
U3 and another half will transform in the H-resonance region to V2m - the second eigenstate 
in matter. Soft decoherence leads to features in supernova neutrinos similar to the ones 
expected within standard oscillations with sin 2 2# 13 < 10~ 4 ; i.e., non-adiabatic conversion 
in the H-resonance [67| . In particular, the earth matter effect should show up. Thus, one 
expects mismatch of the 1-3 mixing measured at the accelerators and bound from studies of 
supernova neutrinos. 

Detailed analysis of the decoherence in matter and effects on supernova neutrinos are 
beyond the scope of this paper. Here we will present a simplified consideration which shows 
that the statistical equilibration is achieved already before the 1-3 resonance. Recall that 
in matter the mass states are not the eigenstates of propagation and therefore oscillate. In 
particular, z/3 should oscillate into a certain combination of v\ and z/2, u a , which depends on 
the density. In our soft decoherence scenario, the decoherence length for supernova neutrinos 
(E ~ (10 — 20) MeV) is rather small: Z/decoh = I/7 ~ 1 m. The oscillation length in matter 
is determined by the refraction length: l m w 2n/V. Therefore for not very large V we 
have Ldecoh "C l m - In this limit, the oscillation effect can be considered in the following 
way. The neutrino trajectory can be divided into intervals of size Ldccoh- A given mass 
state, z/ 3 oscillates on the first interval Z-decoh v z — > a/I — afz/3 + aiz/ a , where |«i| 2 = Pi 
is the transition probability. At the end of this interval the coherence between z/3 and v a 
components of the state is destroyed and in the next interval, L deC oh, they will oscillate 
independently. At the end of the second interval we will have split of the states again, and 
so forth. The probability of transition z/ 3 — > v a in the ith interval, Pj 1 can be estimated 
as 

P % ~ sin 2 2Cass(^) sin 2 fa « sin 2 20£ asB (V i ) 2 « sin 2 29 13 (^L decoh ) 2 , (23) 

where #™ ass = 0^ — 13 is the mixing angle of the mass states in matter (9^ is the mixing 
angle of the flavor states in matter), Vi is the matter potential in the ith interval, fa is the 
half-phase of oscillations in the ith interval . Notice that to derive Eq. (23) we have used the 
constant density approximation for the oscillation probability within each interval of size 
-^decoh- We have then used fa 1 which follows from Ldecoh <C l m . 

Since the initial flux is mainly composed of z/3, the flux transition z/3 — > u a will dominate 
over the opposite transition u a — > z/3, and eventually this will lead to the equilibration of the 
z/3 and v a fluxes. The total transition probability after passing n intervals is given by 

n n 

P « E P < = E sin2 2 ^3^Loh • (24) 

% i 
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This formula is valid when P ^ 1 (linear regime) for which the inverse transition still can 
be neglected. Still the condition P ~ 1 allows to evaluate the length (numbers of intervals) 
over which the equilibration is achieved. Substituting summation in (1241) by integration we 
obtain 

frR 

P^sin 2 29 13 L decoh V 2 (r)dr, (25) 

Jro 

where vr ~ 10 7 m is the radius of the resonance layer. (Recall we are estimating the distance 
from the resonance layer in the direction of center of a star on which equilibration is reached.) 
Taking V = Vji(r^/r) 3 we obtain from fl25|) 

P « sin 2 26 13 [L dccoh V(r )] 2 — ^— . (26) 

decoh 

Taking Vr ~ Am 2 3 /E u , from Eq. fEZS]) . we find that for the whole range of #13 within our 
scenario [see Eq. (|20l)]. propagating from r ~ O.Itr to tr the equilibration condition (P ~ 1) 
is fulfilled. 



V. CONCLUSIONS 

Clearly there is no simple explanation of the LSND result which is consistent with other 
neutrino data. According to some proposals it requires the combination of two exotic mecha- 
nisms. One can ask if any yet unknown mechanism exists which can provide a description of 
(reconcile) all the data. One can take the bottom up approach and try to uncover properties 
(energy and distance dependence) of this unknown mechanism. Our proposal is essentially 
along this line. 

We have proposed an explanation of the LSND signal as manifestation of the quantum 
decoherence of mass states associated to the 1-3 mixing and mass splitting. The decoherence 
leads to a complete or partial damping of the interference terms in the oscillation proba- 
bilities. Our phenomenological scenario is based on the three-neutrino framework (without 
any sterile neutrinos) and makes use of a single decoherence parameter 7 with a sharp en- 
ergy dependence (soft decoherence). The main features of our scenario, which allow us to 
reconcile the LSND signal with the results of other experiments are: 

• zero or negligible decoherence effect on the 1-2 mixing and splitting (712 = 0); 

• decoherence of only the u 3 mass eigenstate; 

• a strong decrease of the decoherence effect with the neutrino energy: 
7 = 7i3 = 723 oc E~ r . 

The strong decrease of 7 with energy allows us to accommodate the LSND signal, while 
being consistent with the null-results of experiments at higher energies, such as MiniBooNE, 
CDHS, NOMAD, and NuTeV. At the same time this energy dependence guarantees standard 
neutrino oscillations in the atmospheric and MINOS experiments. The lower bound on the 
energy exponent, r > 4, follows from the low energy atmospheric neutrino data; r = 4 can be 
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considered as an optimal reference value. The assumption 712 = ensures that oscillations 
in the 1-2 sector relevant for solar neutrinos and the KamLAND experiment are not affected. 

The LSND signal is determined by the mixing angle #13, and the degree of decoherence. 
It implies a lower bound on 1-3 mixing, sin 2 2013 > 0.006 (3er), which corresponds to com- 
plete decoherence. This in turn, leads to testable predictions for upcoming experimental 
searches of 1-3 mixing. For reactor experiments (MeV energies) like Double-Chooz or Daya 
Bay we predict full decoherence already after a few centimeters, leading to no oscillation 
effect when results from near and far detectors are compared. In contrast, the long-baseline 
accelerator experiments (GeV energy range) like T2K or NOz/A are practically not affected 
by decoherence and a signal for #13 should show up. Hence, a mismatch in the #13 mea- 
surements of upcoming reactor and accelerator long baseline experiments would be a clear 
indication for the proposed scenario. At low-energy long-baseline experiments such as the 
CERN SPL superbeam or a Beta Beam with a relativistic 7-factor of 100, our decoherence 
scenario will lead to a distinct energy spectrum of the appearance as well as disappearance 
signals. The soft decoherence can also show up as a distortion in the energy dependence of 
disappearance probability in the second phase T2K. This scenario can also partially account 



for the anomaly found in Gallium radioactive source experiments [55], though explaining 
the full effect might be difficult. Soft decoherence can also affect the supernova neutrinos. 
Despite relatively large #13, decoherence leads to a neutrino composition similar to the case 
of non-adiabatic conversion which takes place in the case of the standard oscillation with 
sin 2 2#i3 < 10~ 4 . Thus, comparing 6*13 measured by T2K and NOz/A and supernova bound 
on this mixing angle can be considered as another way to test the present scenario. 

The energy dependence of the decoherence effect has to be explained by an underlying 
theory that gives rise to quantum decoherence. In this paper we considered the simplest 
power law dependence of 7 in the whole energy range. In general 7 may have a more 
complicated dependence. Indeed, the fact that 7 — ► 00 for E v — > indicates that there 
should be some low energy cut-off below which the dependence of 7 on E v becomes modified. 
The restrictions on 7 come from neutrino data in the energy range from ~ 10 MeV to multi- 
GeV. With a general energy dependence, it is possible that for E u < 10 MeV, 7 remains 
constant or even becomes small again. Such a behavior could modify our predictions for 
low energy experiments, especially reactor neutrinos. Note, however, that under the power 
law assumption the coherence for reactor neutrinos is lost already within a few cm, whereas 
the close detectors are at several hundred meters. So, if 7 is constant below LSND energies 
or even decreases not very fast, we still will have decoherence. Hence, most probably at 
least partial decoherence will be observed and in this case still 1-3 mixing will be different 
in reactor and accelerator experiments. Only in the case of a very sharp cut-off below 
E v ~ 10 MeV (a behavior which appears quite unnatural) one will see the same 1-3 mixings 
in both cases. In any case the decoherence can be still probed by studies of the energy 
spectrum in the phase II of T2K, and CERN beta beam and SPL experiments. Moreover, 
we still expect a disagreement between (9 13 measurements by T2K and NOz/A and the bounds 
from supernova data. 
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